Quantum simulation of spin ordering with nuclear spins in a solid state lattice 
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An experiment demonstrating the quantum simulation of a spin-lattice Hamiltonian is proposed. 
Dipolar interactions between nuclear spins in a solid state lattice can be modulated by rapid radio- 
frequency pulses. In this way, the effective Hamiltonian of the system can be brought to the form 
of an antiferromagnetic Heisenberg model with long range interactions. Using a semiconducting 
material with strong optical properties such as InP, cooling of nuclear spins could be achieved by 
means of optical pumping. An additional cooling stage is provided by adiabatic demagnetization 
in the rotating frame (ADRF) down to a nuclear spin temperature at which we expect a phase 
transition from a paramagnetic to antiferromagnetic phase. This phase transition could be observed 
by probing the magnetic susceptibility of the spin-lattice. Our calculations suggest that employing 
current optical pumping technology, observation of this phase transition is within experimental 
reach. 

PACS numbers: 76.60.-k, 75.10.Dg, 03.67.Lx 
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I. INTRODUCTION 



Many-body problems are common in condensed mat- 
ter physics. However, computations of interesting quan- 
tities, such as critical exponents of a phase transition, 
are difficult, specifically because the dimension of the as- 
sociated Hilbert space increases exponentially with the 
system size. Recent quantum simulation experiments us- 
ing cold atoms in an optical latticeii^ have opened new 
possibilities for the study of these systems. The idea con- 
sists of controlling the interactions inside a many-body 
system so that they take a desired form, and performing 
measurements on the macroscopic behavior of the sys- 
tem. It is interesting to explore how this quantum simu- 
lation concept could be implemented in different physical 
systems. 

Spin- lattice models are a convenient, simple model to 
describe magnetic phenomena^. It was recently realized, 
however, that various apparently unrelated problems can 
be cast into this same language. An example is the 
half-filled Hubbard model in the limit of large positive 
on-site Coulomb repulsion U—. This observation makes 
spin-lattice models interesting not only for fundamental 
studies of magnetism, but also as a means to attack sev- 
eral other problems. Although spin-lattice models have 
been studied for about a century, exact results only exist 
for special cases of low spatial dimensionality. Approxi- 
mation methods such as mean-field theory offer an alter- 
native, but they typically involve assumptions that are 
not rigorously justified. On the other hand, numerical 
methods such as Monte Carlo simulation^ have limited 
efficiency for calculations in real three-dimensional mag- 
netic systems. This motivates the investigation of other 
techniques to study spin-lattice models. 



In this paper, we consider a solid-state lattice of nu- 
clear spins as a simulator of a specific spin-lattice model, 
in which interactions between spins are modulated by 
rapid radio- frequency pulses. The potential of nuclear 
magnetic resonance (NMR) techniques''^ for simulating 
an artificial many-body Hamiltonian is theoretically ex- 
plored. Nuclear spins interact predominantly via a mag- 
netic dipolar interaction, which can be modulated by 
varying the orientation of the crystal relative to the ap- 
plied magnetic field, applying NMR pulse sequences, or 
using different material systems. By optically cooling 
the nuclear spins^i^, or by implementing a similar dy- 
namic polarization technique, it is possible that the sys- 
tem eventually reaches a nuclear spin t emp erature^ i^°i^^ 
where a phase transition-^ occurs. In particular, in Sec. 
|V]we describe the transition from a paramagnetic to an 
antiferromagnetic phase that we expect to occur in our 
proposed system. In this context, we are not consider- 
ing a generic spin-Hamiltonian simulator, but rather a 
problem-specific machine, where we have the ability to 
change some aspects of the spin-lattice Hamiltonian and 
the spin temperature, using the freedom provided by the 
particular experimental system. 

The proposed experiment proceeds as follows: We first 
optically pump a sample of bulk InP in low physical 
temperature and high magnetic field, so that the re- 
sulting nuclear-spin polarization is maximu m^^i^^ . We 
then perform adiabatic demagnetization in the rotating 
frame (ADRF) to transfer the Zeeman order to the dipo- 
lar reservoir-l^iiS,. This is a cooling technique analogous 
to slowly switching off the dc magnetic field in the lab- 
oratory frame. If it is performed adiabatically, i.e. with 
no loss of entropy, the ordering due to alignment of the 
spins along the dc magnetic field is transferred to or- 
dering due to local correlations of the spin orientations. 
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By applying a suitably chosen NMR pulse sequence, we 
can transform the inherent Hamiltonian of the system to 
the desired one in an average sense, to be described be- 
low. In particular, we use NMR techniques to engineer 
an antiferromagnetic Heisenberg model with long range 
interactions. This transformation is validated by aver- 
age Hamiltonian theory (AHT)^^. Our system has two 
species of nuclear magnetic moments. We argue that 
in the ground state of the described Hamiltonian each 
nuclear species is ordered antiferromagnetically, so that 
each nuclear species consists of two sublattices of oppo- 
site spin orientation. Antiferromagnetic ordering with 
two sublattices can be observed in the NMR spectrum 
as the appearance of two peaks of opposite phase: In 
a system with two nuclear species A and B of widely- 
separated Larmor frequencies, we can rotate the spins 
A to the horizontal plane using a ^-pulse, and observe 
their free induction decay (FID) signal. If the spins B are 
in an antiferromagnetically ordered state, they create a 
spatially oscillating local field that modulates the Larmor 
frequency of spins A. This means that spins A in different 
sublattices will produce the FID signals of different fre- 
quencies. However, to decide whether a phase transition 
occurs, we need to measure the magnetic susceptibility 
of the systemilii^. In particular, the longitudinal suscep- 
tibility has a maximum at the transition temperature, 
while the transverse susceptibility exhibits a plateau in- 
side the ordered phase. Nuclear spin susceptibilities can 
be easily extracted from the Fourier-transformed FID sig- 
nals recorded after a short pulse. This measurement ex- 
ploits the Fourier-transform relationship between tran- 
sient and steady-state response of a system to a small 
perturbation. 

All individual pieces of this experiment have already 
been demonstrated, but putting them together is chal- 
lenging. Nuclear ordering has been observed in a series 
of experiment o^^i^^'^°i^^ . However, the effect of NMR 
pulse sequences to engineer an artificial Hamiltonian has 
not been studied in detail, and this is an issue that the 
current study addresses. We should note that the class 
of Hamiltonians that can be simulated with this method 
is severely limited by the characteristics of the particular 
material system, although many interesting cases can be 
studied, for example frustrated spin lattices. 

This paper is organized as follows: In Sec. |TT1 we briefly 
describe adiabatic demagnetization in the rotating frame 
(ADRF), which is essentially a cooling technique for nu- 
clear spins. Sec. |IlT] discusses the calculation of thermo- 
dynamic quantities for systems obeying the spin Hamilto- 
nians of interest, and in particular, the evolution of spin 
temperature during ADRF in our proposed material sys- 
tem. In this way, we can calculate the spin temperature 
reached by the spin- lattice after ADRF. In the subse- 
quent two sections, we describe the spin Hamiltonians 
that can be simulated by this method. In particular, a 
suitable NMR pulse sequence is discussed in Sec. IIVI and 
an experimental material system is presented in Sec. [V] 
Finally, in Sec. IVIl we give an estimate, based on mean- 



field approximation, of the initial polarization needed to 
observe a phase transition. 



II. DESCRIPTION OF ADIABATIC 
DEMAGNETIZATION IN THE ROTATING 
FRAME (ADRF) 

The term adiabatic is used in two different mean- 
ings in physics, one quantum mechanical, and the other 
thermodynamioii^. ADRF is based on the thermody- 
namic definition, which is a reversible process with no 
heat transfer, in which the entropy of the system remains 
constant. The term isentropic also applies to this case. 
By contrast, in quantum mechanics, the term adiabatic 
usually means that the relative populations of the vari- 
ous energy eigenstates of the system are kept constant. 
The two definitions are in general incompatible. 

The system of nuclear spins interacts only weakly with 
the crystal lattice. The relevant timescale for this in- 
teraction is Tr . So, if we are interested in phenomena 
much faster than Ti, we can treat the system as being 
effectively isolated, and describe it in terms of its nuclear 
spin temperatureii, which can be very different from the 
lattice temperature. Of course, if we let the system to 
relax over several timescales Ti, these two temperatures 
eventually become equal. 

Consider a lattice of nuclear spins in a dc Zeeman mag- 
netic field Bq along the z-direction. During ADRF, we 
apply a rotating radio- frequency (rf ) magnetic field in the 
x-y plane of magnitude Bi and angular frequency ui , so 
that the Hamiltonian in the frame rotating with angular 
frequency tui (rotating frame) is 



h[{j,Bo-LUi)i: + j,Bj^ +n'a, (1) 



where I^^y^^- and 7^ are the spin operators and gyromag- 
netic ratio of the i'^ spin dipole, and Ti'jj denotes the 
secular part of the dipolar interaction Hamiltonian. The 
latter follows after dropping the rapidly oscillating terms. 
The three terms in the above equation we called the Zee- 
man, the rf (radio frequency), and the dipolar part of 
the total Hamiltonian, respectively. In the case where 
the lattice consists of nuclear spins of the same species^^ 



t 3 2 



ifi^ - -1+17 I = n'li. 



2 ' ^ 



(2) 



The distance between the spins i and j is r.y, 9ij is the 
angle between the vector and the direction of the 
Zeeman magnetic field z, while if' = if ± z/f and 7/ 
is the common gyromagnetic ratio of the spins^S. Note 
the role of the term — h'fiBiIf in ^ allowing the 
Zeeman and dipolar reservoirs to exchange energy. In 
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the absence of the rf field, the Zeeman and dipolar parts 
of the Hamiltonian are separate constants of motion, and 
the conditions of thermal equilibrium are not satisfied. In 
the presence of the rf field, a common spin-temperature 
is established, and the system is described in the rotating 
frame by the density matrix 



exp 



Tr 



exp 



(3) 



The process of establishment of thermal equilibrium is 
described by the Provotorov equations^^. When the rf 
term becomes comparable to the Zeeman term, Zeeman 
order is established along the direction of the effective 

magnetic field (^Bq — z+Bix. After wi reaches 7/ i?o, 
we switch off the rf field Bi adiabatically, so that all the 
initial Zeeman order is transferred to the dipolar reser- 
voir. 



III. EVOLUTION OF SPIN-TEMPERATURE 
DURING ADRF 

During ADRF, the entropy of the system remains con- 
stant. So, if we have an expression for the entropy as a 
function of the inverse spin-temperature (3, we can calcu- 
late the spin-temperature during the ADRF. In particu- 
lar, we can estimate the desired initial polarization, that 
is required in order to observe a phase transition. 

In the rest of the paper, we drop the superscript rot, 
as we are only going to work in the rotating frame. We 
also use the word temperature to denote nuclear spin- 
temperature. 

In the high temperature approximation-^., we can ex- 
pand the expression for the density matrix ([3]) to first 
order in inverse temperature (3 



P - 



A 



(i - pn) 



(4) 



where A\s a, suitably chosen constant, so that Tr(/5) — 1. 
The entropy S then follows from the equation 



S = —kB^T: (plnp) 



(5) 



For N spin-j particles of the same species, with gyromag- 
netic ratio 7, interacting through the dipolar interaction 
of eq. ^ we find, in the rotating frame 



S 



Nks 
where 



= ln(2j + l)-/32;j2^2j(j±l) 



tin 



Bi 



Bl 
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(6) 
(7) 



In the above equation, we have ignored the rf field Bi. 
Bl can be interpreted as a "local field" due to the neigh- 
boring spins. 



We follow the procedure described in Refs. [2^l23ll2l 
(see also Ref. [2^ for similar calculations) to extend this 
description to the case of lower temperature. The idea 
is to include more terms in the expansion ([4]), applying 
a diagrammatic technique similar to that used in per- 
turbation theory. The Zeeman term is considered as the 
free part of the Hamiltonian, and is treated to all orders. 
The dipolar interaction term is considered as the per- 
turbation, and by grouping the resulting terms using a 
procedure that recalls Wick's theorem^^, we are able to 
efficiently perform the summations. In this calculation, 
we ignore the rf term in ([1]), as this has only the effect 
of bringing the Zeeman and dipolar reservoirs to equilib- 
rium. Also, far from resonance, the two spin species can- 
not exchange energy, because their interaction takes the 
form ^-j, Uifc/f S'f, which cannot induce flip-fiops. When 
the rf part in equation ([T]) becomes comparable with the 
Zeeman part, the two reservoirs can exchange energy, be- 
cause in this case the Zeeman ordering is not along the z 
direction. This means that the Zeeman and interaction 
parts of the Hamiltonian do not commute, so they are 
not separately constants of motion, and the conditions 
for the establishment of thermal equilibrium are satis- 
fied. In our approximation, we do not study the effect of 
different spin-temperatures. Our calculation consists of 
defining the initial polarization of one spin species, and 
calculating the spin-temperature in the rotating frame 
during the ADRF, with the constraints that the system 
is described by one spin-temperature, and that the en- 
tropy is invariant. 



To study the feasibility of observing a phase transition, 
we estimate by the condition that the longitudinal 
magnetic susceptibility vanishes, as the ordered state is 
expected to be an antiferromagnet. This is not an accu- 
rate calculation, however, because only a few diagrams 
are involved in the high-temperature expansion, while in 
the critical region higher order terms become increasingly 
important. A common procedure to take the latter into 
account is described in Ref. 12, in which extrapolation 
procedures are used to compute the limiting behavior of 
high-temperature expansion series given only the knowl- 
edge of the first few terms. Also, antiferromagnets have 
a more rich phenomenology than this simple picture. In 
particular, Refs. HUSHH imply that ID Heisenberg anti- 
ferromagnetic spin chains, in the limit of large spin, are 
gapless for half-integer spin, so that the susceptibility 
does not vanish at the critical point. The mean-field 
descriptionii is quite different, as well. A much bet- 
ter way, based on mean-field theory, to calculate critical 
quantities, is presented in Sec. IVII 
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A. The longitudinal case 

Assume that the spin interaction is longitudinal; that 
is, the Haniiltonian has the form 



1 



u,J'n -Hz+H'. 



2 "y^ i ^3 



(8) 



2J 



A stands for the detuning h (7/ Bo — ^^i)- After thermal- 
ization, the system is described by the density matrix 



exp (^-PoHz - Pn' 



Tr 



exp 



(-/3o 



Hz - f3n' 



(9) 



We allow for different temperatures for the Zeeman and 
dipolar reservoir for convenience in the calculation. In 
the end, we are going to set these temperatures equal. 
Defining the zeroth order thermal average for an operator 
Q as 



(Q>o 



Tr 



exp 



[-PoHz) Q 



Tr 



exp 



{-mz) 



(10) 



and treating the interaction term as a perturbation, we 
have the following expansion for the true thermal average 



(S) = 



E^=o (-/?)" ((^') 2>o/ 



(11) 



^7i'^ involves multiple sums, and the bookkeeping be- 
comes complicated, unless we define the semi-invariants 
M? 



E 



MP 



(12) 



11+12 + ■■■+ip='n 



With this definition, we can write down a diag rammatic 
expansion, like the one described in Ref. [2J. The ad- 
vantage over the brute-force expansion pT|) is that the 
resulting sums are unrestricted, i.e. no constraints such 
as i ^ j are involved in the summations. If we define con- 
nected diagrams as those involving Q, and disconnected 
ones as those not involving it, we can factor the numer- 
ator as the product of the sum of connected diagrams 
times the sum of disconnected ones. The latter serves to 
cancel the denominator, so that only connected diagrams 
matter. 

The calculation of the semi-invariants is straightfor- 
ward. With the definitions 



V 



/,(77)=tr(e''^") 



_ sinh [vjj + ^)] 



sinh 



(13) 



the first few semi-invariants can be written as 



,(2) 



= t 



3 

(3) 



(4") . 



(fT 



(14a) 
(14b) 

(14c) 



The general case is covered in Appendix [X] 



B. Calculation of thermodynamic quantities 

In this subsection, we calculate the entropy as a func- 
tion of inverse temperatures P and /3o and detuning A. In 
the final formulas, we always set Pq = (3, but we consider 
P and Po as independent variables in the calculation of 
derivatives. We first need to calculate the average polar- 
ization. 



dp^ 



_d_ 

Wv. 



d 



{-H'l 



drj 



(15) 

The dipolar interaction energy (7i^) is calculated in Ap- 
pendix |B] in terms of the functions Gn(.) defined there. 
Performing the integral in the above equation, we find 



{h{Po,p) ^ ^(1) 



N 

E 



p 



n+l 



d 



n + 1 dij 



G 



„(t«(r;),t(2)(r;),...,i("+i)(^) 



(16) 



If we keep terms to first order in A, we can estimate the 
longitudinal magnetic susceptibility. 



1 

TV 



-iPA) 



jU + 1) 

540 



180 - (33 + 56j + 56f)P^ {jhBLY 



(17) 



This formula provides a method to estimate the criti- 
cal temperature of the paramagnetic to antiferromagnetic 
phase transition from the temperature at which the sus- 
ceptibility becomes zero. For spin-^ and spin-|, we find 



ksTcU = -) = 0.645 {jHEl) , 



ksTcij 



2.81 (jhBL) . 



(18a) 



(18b) 



This is a very rough estimate, as discussed in the be- 
ginning of this section. Further, it refers to the ordering 
of each species in its own lattice, not involving coupling 
between them. Eventually we will use these couplings to 
create the artificial Hamiltonian described in Sec. IIVI 
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The free energy follows from the expression 

lnZ(/3o,/3) =lnZ(/3o,0)- / {n'D){po, X)dX, (19) 

Jo 

which can be expressed in terms of the functions G„ 

^lnZ(/3o,/3)=ln/j(-/3oc^o) 
- E (t(i)(r;),t(^)(r,), . . . , . (20) 

n 

We are now ready to calculate the entropy. 

^ = PoL^oih + Pin'o) +lnZif3o./3); (21) 
Kb 

in terms of the functions Gn-, 

^ ^ -,4" + 1" /, W + S ^ + "|c,.) ^ 

(22) 

C. Results 

During ADRF, the entropy remains constant. In order 
to obtain an ordered final state, we need the initial state 
to be sufficiently polarized. A rough estimate for spin-^ 
species is to set the spin temperature ksT of the ini- 
tial paramagnetic state equal to the level spacing h^BQ. 
This means that we need the initial polarization, which is 
defined as the thermal average magnetization over max- 
imum magnetization, to be 

Sec. IVII gives a better picture of how the initial polariza- 
tion influences the ordering observed under the artificial 
Hamiltonian we wish to implement. 

In Fig. [1] we present the results of our numerical cal- 
culation. About 50% initial polarization of the spins is 
required in order to reach a spin temperature where a 
phase transition can be observed. Details about the pro- 
posed material system are presented in Sec. fVland in the 
Appendices. 



IV. PULSE SEQUENCE AND THE ARTIFICIAL 
HAMILTONIAN 

If the system has been sufficiently cooled via opti- 
cal pumping and ADRF, it might be in dipolar-ordered 
state because of the inherent interaction between spins. 
On the other hand, by applying a WAHUH A^^i^" , or 
MREV8^i^i2^ pulse sequencei^ on resonance, we can 
decouple spins of the same species, and transform the 



Evolution of temperature during demagnetization 
spin j = 1/2, initial polarization = 40% 
3.5| ' ' • • 
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FIG. 1: Evolution of spin temperature T vs detuning A during 
ADRF. Dotted line: high temperature approximation. Bro- 
ken line: first order approximation. Continuous line: second 
order approximation. These approximations are discussed in 
Sec. mil and in the Appendices. Insets: Initial polarization of 
50%. The parameter Bl (local field) is defined in (O, Tc is 
calculated in p8|l . The calculation involves one spin species 
with spin-i or | in an fee lattice and the Zeeman magnetic 
field along the (111) crystalline direction. For high initial 
polarization, large spin, and near resonance, the calculation 
does not converge. 

interaction between unlike spins from an Ising-like inter- 
action to a Heiscnberg-like interaction. 

In particular, the WAHUHA pulse cycle is given by 

T-P,-T-P^y-2T-Py-T-P^,- r, (24) 

where P^ and P-a denote ^ and — | pulses respectively 
around the a-axis, (a = x^y). t and 2t denote the time 
interval between the application of consecutive pulses. 
The pulse cycle time is tc — 6t plus the duration of 
the four pulses. MREV8 contains two WAHUHA-type 
pulse cycles, concatenated to mitigate pulse imperfec- 
tion effects. If the pulses are short enough, then their 
effect can be modeled as an instantaneous transforma- 
tion of the Hamiltonian. To first order in tc relative to 
the magnitude of the Hamiltonian, the system dynam- 
ics are described by the time-averaged Hamiltonian over 
the pulse cycle. Average Hamiltonian theory (AHT) is 
covered in detail in several textbooks (see, for example, 
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Ref.[ll). To first order in AHT, the dipolar interaction 
in the rotating frame between nuclei of the same species 
([2]) is eliminated. If we simultaneously apply this pulse 
sequence to both species, their dipolar interaction in the 
rotating frame 



i,k 

is transformed to a Heisenberg-like interaction 



T-i'is = ^Wik^ii ■ Si 



(26) 



i,k 



In the above sums, the index i runs over the spins /, 
whereas the index k runs over the spins 5*. Wik has a 
similar expression as in 

Ha 1 - 3cos^(6l,fc) 
Wik = -;-n ^ns 5 • (27) 

rfk 

If we use the modified WAHUHA pulse cycle 

ri-P^-T2- P-y - 2t2 -Py-T2- P-^ - Ti, (28) 

the effective interaction Hamiltonians to first order in 
AHT are 



H'/7(ri,r2) 



Tl - T2 



n'ls{Tl,T2) 



i,k 



Wik 



2t2 

t 



— {lfS\ 



I'. SI 



(30) 

where = 2ri +4r2. By keeping tc constant and varying 
T2 adiabatically from to -^^ , we can get an interpolation 
between the initial and final Hamiltonians 7i™'*, 
and transform the state of the system from the thermal 
state of Ti™^* with entropy S to the thermal state of Ti-^*" 
with the same entropy. In this way, we can study the 
thermodynamics of the particular Heisenberg model for 
varying entropy. An alternative method is to keep ri = T2 
and vary the width of the pulses, so that they induce 
rotations by an angle 0, same for all pulses. By varying 
9 from to ^, we have a similar effect, but the equations 
for the effective interactions are more complicated. 



quadrupole term in the Hamiltonian, which gives a non- 
zero contribution in the presence of electric field gradi- 
ents. When this contribution is the dominant one, it can 
be expected that any ordering due to the dipolar interac- 
tion term is destroyed. Phosphorous has only one isotope 
•^^P with spin-i, while Indium has two isotopes, ^^^In 
(25) and ^^^In have the same spin (|) and are close in Lar- 
mor frequency for typical NMR magnetic fields, so we can 
consider Indium to be homonuclear. InP is a direct-gap 
semiconductor with strong optical properties, and effi- 
cient optical cooling has already been demonstrate d^^'^'^ . 
InP has a zincblende structure; if a dc magnetic field is 
applied along the (111) crystalline direction, and cou- 
plings between spins of the same species are neglected 
during the application of the NMR pulse sequence as de- 
scribed in Sec. IIVI then the magnetic dipolar interac- 
tion between nuclei creates the following structure (note 
that the gyromagnetic ratio 7 of both In and P nuclei 
is positive) . The planes of Indium and Phosphorous nu- 
clei are grouped into pairs. The coupling between near- 
est neighbors (nn), which are always of different species, 
is antiferromagnetic if both nuclei belong to the same 
pair of planes, and ferromagnetic if the nuclei belong 
to successive pairs of planes. Second- nearest-neighbors 
are of the same spin species, and so their interaction is 
averaged out. Some third-nearest-neighbors interactions 
cause frustration, but the magnetic field orientation en- 
sures that their effect is small. Ignoring interactions far- 
ther than third-nearest-neighbors, this leads the rather 
nice picture shown in Fig. [21 Interactions within each 
pair of planes are antiferromagnetic. Coupling between 
plane pairs is ferromagnetic. It is possible (but by no 
means certain) that in the absence of an offset dc mag- 
netic field in the rotating frame, each plane pair acts as 
a two-dimensional antiferromagnet, where the ferromag- 
netic interaction between plane pairs merely ensures that 
all of them have the same structure. A complementary 
picture results if we consider each spin species indepen- 
dently. In our proposed structure, each species forms a 
lattice with antiferromagnetically ordered spins. Anti- 
ferromagnetic ordering within spins of the same species 
results in an indirect way through interactions with spins 
of the other species, as shown in Fig. [2l 



How reasonable is it for us to neglect interactions be- 
yond 3rd nn? The strength of dipolar interactions falls 
off with distance only as p-, and because the volume in- 
tegral 



V. CHOICE OF THE EXPERIMENTAL 
SYSTEM 

Bulk InP is a suitable material for this experiment. As 
described in Sec. IIVI two spin species are required to 
implement the desired model Hamiltonian. A strained 
lattice, on the other hand, like a quantum well ex- 
hibits strong quadrupolar broadening. This is because 
spins higher than i, such as ^^^In and ^^^In, include a 




is logarithmically divergent, interaction with distant nu- 
clei, could significantly modify the above picture. We 
discuss this question in the next section, in which we 
take into account coupling between every pair of spins 
within the mean-field approximation. 



7 



spin as being subject to a local field 




-2 



FIG. 2: (color online) The configuration of the nuclear spins 
under the experimental conditions described in Sec. |V] 
Spheres of the same color represent nuclei of the same species. 
Solid lines are antiferromagnetic bonds, while broken lines are 
ferromagnetic ones. Only nearest-neighbor bonds are plotted. 
The perpendicular direction in this figure is the (111) crys- 
talline direction. We can see the postulated antiferromag- 
netic planes, which are coupled ferromagnetically. Within 
each plane, spins of one species tend to form a ferromagnetic 
lattice which is aligned antiferromagnetically to the lattice 
formed by the spins of the other species. On the other hand, 
spins of the same species belonging to successive planes tend 
to align antiferromagnetically. 



VI. MEAN FIELD THEORY OF ORDERING 

In Sec. mil we presented an estimate of the initial nu- 
clear spin polarization needed in order to observe a phase 
transition from a paramagnetic to an antiferromagnetic 
phase. Our discussion involved thermodynamic argu- 
ments. A more exact estimate could not be obtained 
with this method, as our calculation includes only a few 
diagrams, and is not valid near the phase transition. Ad- 
ditionally,It does not converge to an answer when the spin 
temperature is small, and the calculation of Tc is not ac- 
curate, as already explained in the beginning of that Sec. 
mil We now use the local Weiss field approximation to 
study the problem of ordering of nuclear spins under the 
effective interaction of our model. Our discussion closely 
follows Refs. [iHl. 

During the application of the decoupling pulse se- 
quence, the spin Hamiltonian in the rotating frame is 



n = C7/75 



(32) 



where 



l-3cos''(eifc) 

range correlations between spins, we can consider every 



Neglecting short- 



k 



(33) 
(34) 



We then assume thermodynamic equilibrium, so that the 
expectation values (7^), and (Sk) equal their thermal av- 
erages. Near a phase transition, the magnetization is 
small, so that we can neglect all terms except the linear 
term 



(Sk) =--S{S + l)/3hLJk. 

(35) 

/ and S are the spin quantum numbers for spins / and 
S. We then write a self-consistent set of equations, using 
the lattice Fourier transforms of (7^), (Sk), and (Aik): 

" A//"(fe) + cjnsA*{k)S"ik) - , . 

c777sA(fc)/"(fe) -I- XsS"{k) = 



+ l)/3, A 



-1 
s 



(3S{S + l)/3, a e 



where Xj 

{x,y,z}, and I°'{k), S°'{k), and A{k) are the lattice 
Fourier transforms of (/f ), (S^), and {Aik) respectively. 

We now calculate the transition temperature and or- 
dered structure. The above equations are satisfied for 
all k at the phase transition. Thermodynamic stability — 
implies that we should set all I°'{k), and S"{k) equal to 
zero, except the one which gives the maximum \Tc\. We 
have to keep in mind that this can occur for small |A;|, in 
which case the predicted ordered structure is ferromag- 
netic with domains. 

For a zincblende lattice, with coupling only between 
different species and the dc magnetic field parallel to 
the (111) crystalline direction, this maximum occurs for 
ko = a-i(f , f , f ); \A{ko)\ = T.OOa-^, with 2a the edge 
of the conventional cubic cell. The predicted ordered 
structure is the same as the one described in Sec. |Vl 
which followed from nearest neighbor interactions alone. 
Because of the Heisenberg form of the interaction, there 
is complete rotational degeneracy for classical spins. The 
transition temperature is calculated to be 



knTr 



c7/7s|^(feo)|- (37) 



The spin temperature of the system is modified by the 
ADRF process, and the subsequent application of the 
pulse sequence. So, the relevant quantity is the critical 
entropy, or the critical polarization in high field. To cal- 
culate it, we follow the formalism described in Sec. IIIII 
To avoid complications, we use the high temperature ap- 
proximation. For a lattice with N spins / and N spins 
S, the entropy before ADRF is 



= ln(2/+l)-^/3|7|ft2^^)-^^ 



+ ln(2S'+l) 



1 

2' 

^/3|7i^^^^(i3o)^(38) 
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while the polarization of spin-i nuclei (^^P for InP) is 

p.tanh(^). (39) 

We first assume that both spin species always have the 
same spin temperature. This is not true for the opti- 
cal pumping mechanism for species with different spins 
and/or different gyromagnetic ratios, and we study this 
effect later. During the application of the pulse sequence 
on resonance to both species, the entropy is 

^=ln(2/+l)+ln(25+l) 

-^/3^^^^^^(3c^7l7|)E(^-)'- (40) 

k 

The above equation is valid only in the paramagnetic 
phase. For the zincblende structure with the applied Zee- 
man magnetic field Bq \\ (111), Efc(^ifc)^ = 13.238a-6. 
At P = Pc, this expression for the entropy is equal to the 
value in eq. ((38| . If both species are spin-i, and have 
identical gyromagnetic ratios, the critical polarization is 
Pc = 56%. 

The general case does not include any further com- 
plications. Just after ideal optical pumping of the 
sample^i^, the spin temperatures of the two species sat- 
isfy 

Pni = /3s7s- (41) 
Equilibration between the two spin species occurs 
for Bi ^ (^0-^), Bs ^ (^0-^), and 

II [Bo - ^) ~ 7s (Bo - where Bj, Bs are the 

rf fields acting on the / and S spins. The first condition 
arises as we need the interaction term to not commute 
with the rest of the Hamiltonian for / and S spins (see 
ID)), so that the two spin species can exchange energy. 
The second condition maintains that simultaneous spin- 
flips conserve energy, so that they have a non-negligible 
probability. We assume that these fields are much larger 
than the local fields of Eq. (O, so that wc can ignore 
the latter. We also assume that equilibration happens 
at some particular point of ADRF, when the / and S 
spins have temperatures f3j and respectively. From 
equation ([6]) 

P'i(Bo-^)V2 = PjB,^^^P^^^ 
P's[Bo-^)V2^PsBo) P's Psls 

That is, the two species are already in equilibrium, and 
there is no entropy loss. 

It is now straightforward to calculate the critical po- 
larization for both species, with the assumption of ideal 
optical pumping and no entropy loss during the transfor- 
mation of the Hamiltonian. The result is 

pI = 15%, pi" = 49%. (43) 



Conclusion 

We have proposed an experiment of quantum simula- 
tion of spin Hamiltonians by the manipulation of nuclear 
spins in a solid. Bulk InP seems to be a suitable material 
for this experiment. A particular orientation of the crys- 
tal produces a convenient model Hamiltonian. We have 
presented an NMR pulse sequence, which transforms the 
natural dipolar spin Hamiltonian to a Heisenberg model 
with long range interactions. We have also performed a 
calculation to estimate the experimental feasibility of our 
proposal. Our results suggest that initial polarization of 
Phosphorous nuclei of about 15%, and of Indium nuclei 
of about 50% is needed in InP, in order for a phase tran- 
sition to be observed. The difference in these values is 
because of different nuclear spin quantum numbers, and 
this values should be consistent under ideal optical pump- 
ing conditions. This critical polarization seems reason- 
able for optical pumping of bulk InP in low temperature 
and high magnetic field. 

This scheme provides an alternative for quantum sim- 
ulation experiments which is predicted to be technically 
simpler than the cold atom approach, and involves a 
macroscopic number of quantum particles. It is true 
that its flexibility in implementing a specific spin-lattice 
Hamiltonian is limited, and the method of detection de- 
scribed in the introduction cannot give any information 
about the specific orientation of the spin dipoles. To 
extract this information, one needs more complicated 
techniques, such as neutron scatteringi^. More sophisti- 
cated NMR measurement strategies might also be pos- 
sible. Nevertheless, this scheme enables the study of 
several interesting spin-lattice models. Similar experi- 
ments have already been performe d^^'^^i^^ , and they con- 
stitute a proof of principle for the validity of the spin- 
temperature concept, and for this method of manipulat- 
ing the dipolar interaction. However, in these experi- 
ments only one spin species was used, which limits the 
freedom to manipulate the internal Hamiltonian. In par- 
ticular, -l-pulses could not be used, as the interaction 
Hamiltonian is averaged out, and the equilibration time 
becomes infinite. In our case, the interaction between 
different spin species remains non-negligible, and the adi- 
abatic criterion can be satisfied near resonance, when the 
two spin species can exchange energy. 

Our proposal is known to have several uncertainties 
and approximations, discussed below. To our knowl- 
edge, nuclear ordering has not been observed until now 
using optical pumping as the cooling technique. This 
renders the experiment both interesting and technically 
challenging. We also have to point out that this exper- 
iment, like most NMR experiments, involves subtleties 
because of various relevant timescales, and the inequal- 
ities that they have to satisfy. A detailed discussion is 
quite technical, but plenty of information can be found in 
the literature, in particular about when AHT is validly, 
and when ADRF^'* is adiabatic. What is not adequately 
studied, however, and requires further research, is the 
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rate of establishment of thermal equilibrium in low spin 
temperature while a multipulse sequence is applied. This 
scheme also suffers from several non-idealities, namely in- 
direct nuclear magnetic interactions"^^-^^, non-idealities 
of the multipulse sequence, or not perfectly adiabatic 
transformation of the Hamiltonian. Nonetheless, non- 
ideal ADRF is found in Ref. ^3§ not to pose a serious 
problem for high initial polarization and typical NMR 
applied magnetic fields. Further, the method described 
for the calculation of the transition temperature is valid 
only for classical spins within the mean-field approach, 
so a more sophisticated calculation would be interesting. 



The functions t\p' are defined in (fT3|) . 

APPENDIX B: CALCULATION OF THE 
DIPOLAR INTERACTION ENERGY 



As is clear in Sec. IIII Bl calculation of the dipolar in- 
teraction energy is the key to the calculation of other 
interesting thermodynamic quantities, and in particular 
the entropy, which is our ultimate goal. The dipolar en- 
ergy is 
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APPENDIX A: THE SEMI-INVARIANTS FOR 
THE GENERAL DIPOLAR INTERACTION CASE 



(Bl) 

with Vij = — ^. N is the total number of spins. Because 
of the symmetry of the resulting diagrams, the terms 
{I^I^) and {I^I^) give the same contribution, so we 
keep just the former, multiplied by vu. The second order 
expansions give the diagrams of Fig. [3] and result in the 
following expressions 



In Sec. IIII Al we calculated the semi-invariants in the 
case that the spin interaction is longitudinal. We now 
restore the transverse term of the Hamiltonian 



2 *>i 



''villi 



For the dipolar interaction, Vi 



Hz + 'W'd • 
(Al) 

-\uij. We take into 
account that the operators , if do not commute and 
introduce time-ordering. Details can be found in Refs.l2^. 
[2^. The semi-invariants are defined in a similar way as in 
(fT^ . but the partition sets of the right part must follow 
the ordering of the left hand side. The additional semi- 
invariants we need for our calculation are: 



b) 



FIG. 3: Diagrams for the calculation of the expansions in 
equations a) (|B2I) and b) (IB3|) . The diagrammatic rules are 
presented in Refs. HIIH. 



MO (/+/-) 


= .j(j + i) + tr 


y-(2) 
''J ' 




(A2a) 


M°(/-/+) 


= jij + 1) - tf^ 


V-(2) 

' 




(A2b) 


M^{F-i+r) 


= M^ii+rp) -- 










+ f)(i + if))- 


y.(3) 




(A2c) 


Ml{rri+) 


= M^{ri+r-) -- 










-f>(i-,5")- 


.(3) 

' 




(A2d) 


M°{i+rr) 


+ M^M^{i+r) 








+ 1) 




+ 2f) 


^(3) 


(A2e) 


M°{rF-i+) 


+ M^M^{ri+) 








jU + 1) 




-2f) 


^(3) 


(A2f) 



{ill!) {-l3){M'2fuu + {-PY-{Mlf{uuf + 



-(-/3)2MO(+-z)AfO(- 
-(-/3)2(mO)3E^.i,^.,,, 



^)( 



Vli.2 



(B2) 
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(/+/r) ^ (_/3)^M°(+-)Af2"(-+) + {-pffuu^[Ml{+z-)Ml{- 



(B3) 



The final result can be expressed in terms of the following 
quantities, which depend on the particular lattice and the 
orientation of the magnetic field: 



i 

i 



(B4a) 
(B4b) 
(B4c) 



^,J 



For the case of an fee lattice, with the magnetic field 
parallel to the (111) crystalline direction, we obtain the 
values 



h = 1.6908 



\4:TT 



(B5a) 



h = -0.0072073 



= 2.1173 



47r 



(B5b) 
(B5c) 



where 2a is the edge of the conventional cubic cell. We 
refer to these values in Sec. IIIICI 

With these definitions, we have the following expansion 
for the dipolar interaction energy: 



l(?i;i,)=/3Gi(tW(ry),<(^)(r,))+ 



(B6) 



for Vij = —^Uij, the functions Gi, and G2 take the form 



Gi (iW(ry),t(2)(r,) 



(M^fh + -M^{+-)M^i-+)l2 



(B7a) 



liM^'h + z)Mli- + z)h + {Mi^fK, + 



(Af3"(+z-)Af3"(-z+) + A4{z + -)Ml{z - +)) /3 - 1^mO(+-)m"^(-+) {m",{+-) + AI^i-+)) K,] . (B7b) 
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